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We study the entanglement dynamics in the system of coupled quantum fields. We prove that if 
the coupling is linear, that is if the total Hamiltonian is a quadratic form of the field operators, the 
entanglement can only be transferred between the fields. We show that entanglement is produced 
in the model of the two-mode self-interacting boson field with the characteristic Gaussian decay of 
coherence in the limit of high number of particles. The interesting feature of this system is that 
the particles in different modes become entangled even if there is no direct interaction between the 
modes. We apply these results for analysis of the entanglement dynamics in the two-mode Jaynes- 
Cummings model in the limit with high number of photons. While the photon-atom interaction is 
assumed to conserve helicity the photons with different polarizations still get entangled with the 
characteristic entanglement time linearly increasing with the number of photons. 

PACS numbers: 03.67.Bg,03.65.Yz,32.90.+a 

I. INTRODUCTION 

Entanglement is a feature, which is quintessentially quantum. It signifies that several particles may form a new 
entity, a complex of particles, say, not two particles but a pair and so on. If, for example, one would attempt 
to separate the particles from such complex by performing a measurement on one individual particle, one would 
unavoidably modify the states of others even if the direct interaction between the particles is absent or negligible. 
Such departure from the classical properties not surprisingly attracts the significant attention (see e.g. the recent 
extensive review Ref. [ll). 

As the problem of special interest stands out the problem of sources of entangled states. For instance, nowadays the 
most developed and widely used method of generating the entangled photons is the parametric down conversion;^ 
which is based on on the two-photon radiative decay of excited states of a matter. This method, however, suffers from 
intrinsic limitations — very low yield and rescaling the wavelength of the emitted photons4r— Therefore, there is the 
constant search of alternative sources of entangled light. 

The consideration of the problem of entanglement of photons in the course of interaction with the matter excitations 
naturally suggests the more general perspective of the systems with non-conserving number of particles, which are 
described in the field theoretical framework. In this approach the particles appear not as a predefined entities as, 
say, qubits in canonical quantum mechanical treatment, but rather as the excitations of the quantum field. In this 
context the difhculty of producing the entangled states is clearly demonstrated by the simple example of a bosonic 
field excited by an external source. Let the field dynamics be described by the Hamiltonian 



H 



XI (efe4afc + efc(04 + 4Wafc) , (1) 



where k enumerates the modes of the field with the energies e^, a\. and are the creation and the annihilation 
operators of the particles in the respective modes and are c-numbers determined by the projection of, generally 
speaking time dependent, external classical field on the modes. In order to avoid complications arising in the fermion 
casei""— in the present paper we deal with entanglement of only boson fields. 

Before applying rigorous methods for analysis of this situation let us first consider it using often employed arguments 
based on the interference of different paths leading to the final state of the field. Let us assume for the sake of this 
discussion that either the system admits only two modes or for some reason we can restrict ourselves to considering 
only two of them out of the whole set. Thus the sum in Eq. ([T]) would contain only terms with k = 1,2. The typical 
argument sounds as follows. One path is when one particle first goes under the action of the external source to mode 
1 then to mode 2, while the sequence of the modes for the second particle is — !■ 2 — 1, where denotes vacuum. 
The second path for the pair is obtained by "flipping" the particles: — 2 — 1 for the first particle and — )■ 1 — )■ 2 
for the second one. The interference between these two paths is assumed to lead to entanglement. 

More rigorous description of entanglement uses the notion of the single-particle observables. Each such observable 
is represented by a single-particle operator O = Okqa\.aq (see e.g. Ref. [l^), which motivates introducing the 
single-particle density matrix (SPDM) 



kqit) - (4(t)a,(t)) . (2) 
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Throughout the paper we incorporate the time dependence into the Heisenber representation of the field operators 
flfe(0 ~ exp(i'Hi)afc exp(— zHi), therefore, the average in formulas similar to Eq. ^ is taken with respect to the initial 
state (...) — {ilj{t = 0)1 . . . — 0)). For the purpose of this consideration it suffices to assume that entanglement 
manifests itself in incoherent single-particle states, that is the SPDM has rank more than one;^*^, The solution of the 
operator equation of motions can be easily seen to have the form 

Ofc (i ) = flfc (0) exp i-iekt) + Gk (t) , (3) 

where Gk{t) = /J dt' exp[— iej.(i — t')]Cki^') ^^'^ c-numbcr functions. Substituting this representation into Eq. © and 
assuming that the initial field state is the cigenstate of the operator of the total number of particles, Af — J2k '^k'^k, 
we express the time dependence of the SPDM in terms of its initial value Kkq{0) and the contribution due to the 
excitation Ekqit) = Gl{t)Gq{t) 

Kkq{t) = e"(^''-^')iffc,(0) + Ekqit). (4) 

Generally the dynamics of entanglement in the presence of the external source is not trivial. The addition Ekq{t) may 
lead to the change of entanglement depending on the structure of the initial state. If, however, initially the system 
is in vacuum state then matrix E{t) (throughout the paper the hats denote matrices in the space of the modes of 
the quantum field) with the matrix elements Ekq{t) is at most of rank one. Hence, if K{Q) — then the rank of the 
SPDM is no greater than one implying the absence of entanglement. 

This simple example shows that the problem of the dynamics of entanglement should be treated with certain care. 
First, the naive arguments based on the picture of interference of different paths connecting initial and final many- 
particle states may be misleading, when it is applied to the field dynamics. One has to study the actual structure of 
the field state. Second, no matter how complex is the spectrum, efc, of the (boson) quantum field the states, which are 
reached out of vacuum under the action of the classical external excitation, are disentangled. Thus, in the discussion 
above regarding the effect of the excitation on an initial entangled state that state should be created by other means 
than the external excitation. 

The rest of the paper is organized as follows. In Section |lT] we consider the entanglement dynamics for the system 
of linearly coupled boson fields. In Section ??e study the entanglement for the two-mode self-interacting boson 
field. In Section IIVI we apply the results obtained for the analysis of the entanglement dynamics in the two-mode 
Jaynes-Cummings model in the limit with high number of photons. 



II. ENTANGLEMENT TRANSFER BETWEEN COUPLED FIELDS 



We start the analysis of the entanglement dynamics from the consideration of two coupled fields described by the 
Hamiltonian 

^ = E 4"^4«fc + E ^« + E (dk,.alb. + dlAak) , (5) 

k K k,K 

where k and n enumerate the modes of the fields, e{,"\ e^^^ and d^.K are the spectra of the fields and the coupling 
constants between them, respectively, and the operators ak and 6^ are assumed to obey the boson commutation 
relations. 

The important feature of the time evolution of entanglement in this system is that the total entanglement remains 
constant and is solely determined by the initial state. More precisely, instead of the operators at and let us 
introduce the combined operators Ui. That is instead of two fields with two sets of modes and Mb we consider 
the single field, whose modes are the direct sum Ma ® M^. Thus Ui = if i G Ma and Ui = 6^^ if i G Mb. For field 
u we define the SPDM 

K,,{t) ^ (u\{t)u,{t)) . (6) 

When both i and j belong to, say. Ma the respective matrix elements Kij{t) give the SPDM for field a and so on. 
Furthermore, let the matrix K have the spectral representation 



(7) 
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where A[ and V[ are the eigenvalues and the unit eigenvectors of K, respectively, and (S) denotes the tensor product, 



which is defined as (v (g) v' 



Then the total entanglement can be found as the von Neumann entropy 



(8) 



where A[(i) = Xi{t)/ J2m ^m{t)- The total entanglement can be shown to remain constant, so that it is determined by 
the initial state -EAr(t) — i?jv(0), and the entanglement evolution restricts purely to its redistribution. 

We prove this in a bit more general context. Namely, we consider the system described by the Hamiltonian 



(9) 



where the field operators obey the commutation relation 



[Uj,Ui\ = Vij 



(10) 



with Vij being a symmetric c-number matrix. The total entanglement is defined as the von Neumann entropy of the 
SPDM. Its conservation follows from the fact that the evolution SPDM is governed by the equation 



L,K 



(11) 



where Lij = hikVkj- Indeed, the spectrum of the matrices, whose time dependence is determined by an equation 
of motion with the commutator in the right-hand side, remains constant. Hence, in the spectral representation similar 
to presented in Eq. ([7|) only the eigenvectors will be functions of time. 

For example, if initially there was only single non-zero eigenvalue (i.e. entanglement was zero) it remains the only 
one later on implying no production of entanglement. The same result holds for the SPDM corresponding to either 
fields a or h. Indeed, the SPDM of the field a is obtained from the SPDM of the combined field u applying the 
respective projection operators. Ha, so that 



K 



(a) 



(12) 



Using for K its spectral representation one can see that such projection cannot increase the rank of the SPDM. 

At the same time the circumstance that the time evolution of entanglement of specific particles is determined by 
the projections of the total SPDM results in non-trivial time evolution of initially entangled state. As an example let 
us consider the situation of small total entanglement, more specifically, when there are only terms with [ = 1 and 2 
in Eq. © with Ai > A2. According to Eq. the SPDM for the field a is given by 



(13) 



where v 



(a) 
t 



HaV,/ 



n„v. 



and \[''\t) = A[ 



n„v. 



Thus, the value of entanglement depends on the magnitude of 



< 1 



the projections IlaVi, which is determined by the internal dynamics of the coupled fields. In particular, if 

one may have A^"'' {t) ^ A2"'' {t) resulting in the significant entanglement of particles a, either with each other or with 
particles h. For more concrete information one has to take into account that, generally speaking, Eq. (IT^ may not be 
the spectral representation of the matrix K^""^ because the vectors v^"-* {t) and v^"-* {t) are not necessarily orthogonal. 
Let us denote 9(t) = \^f'f\t) ■ ^^\t)\'^, then the eigenvalues of K^°'^ are found to be 



(a) _|_ y(a) 



1 

± - 

2 



(a) 



+ 4Al"^A^"^0(t). 



(14) 



Thus, if at some particular instant one has \^i\t) ~ then, depending on the details of the dynamics of the 



eigenvectors, one may have E 



(a) 
N 



1. Strong entanglement is produced when the "weak" component of the SPDM 



^ A2 is more effectively transferred into the field a than the major component ^ Ai, whose only small part is moved 
into the field a during the evolution. At the same time, as can be seen from Eq. ([T4|) after slight change of notations, 
such spike of entanglement between particles a is accompanied with disentanglement of particles b. 

The main condition for this geometric effect is the smallness of the projection of the respective eigenvector of 
the total SPDM. As a result the characteristic feature of the SPDM of the strongly entangled states in this case is 



4 



Tr[^(")] <C Tr[^] = Tr[^('^)] + Tr[^(*')]. That is the states with E^j^^ « 1 developed from the states with low total 
entanglement are characterized by low excitation, while the strongly excited field, say the field b in the considered 
example, for which one has Tr[i?(^^] sa Tr[i^], remain only weakly entangled. 

The fact that the entanglement of the specific particles is determined by the projections of the total SPDM may 
lead not only to increased entanglement of the particles but also to disentanglement. Indeed, if for the SPDM given 
by Eq. (O with I = 1,2 at some instant the vectors vi and V2 belong to different subspaces (say, IlaVi = vi and 
nf,V2 = V2) then the entanglements of both particles a and b are zero. One can understand this effect introducing the 
isospin quantum number, so that one value of isospin corresponds to the particle a and another one stands for the 
particle b. Such situation, when non-zero total entanglement coexists with zero entanglement as seen from specific 
particles' SPDM, corresponds to the entanglement stored in isospin. 

It should be noted that while we have limited our attention only to boson fields the same formal results can be 
easily checked to hold also for fermion, i.e. obeying the anticommutation relations, fields. 



III. ENTANGLEMENT PRODUCED BY INTERACTION 



The consideration in the previous section showed that simple linear coupling of the quantum fields is not enough 
for the fields to become entanglement. The situation, as we will show in the present section, is different when there is 
an interaction in the system, that is when the energy depends non-linearly on the number particles. More specifically, 
we consider the self-interacting two-mode boson field described by the Hamiltonian 

"H = ^ Efca^afc + ^ Ukalalakak + U+-a''_al_a+a^. (15) 
k=+,- fe=+,- 

Here the first two terms describe the internal dynamics of the modes with ek and Uk being the energies of the modes 
and the interaction parameters, respectively, and the last term represents the interaction between the modes, which 
is chosen in the form admitting conservation of particles within each mode. 

In order to evaluate the entanglement we need to inspect closer its relation to the matrix elements of the SPDM 
and the state of the field. We start from the case of a two-particle case, which provides the clear connection with the 
standard quantum-mechanical consideration. In the basis of the population numbers any two-particle state can be 
presented as 

\^) = a2fi |2, 0) + aia |1, 1) + ao,2 |0, 2) , (16) 

where denotes the state with n+ particles in the "-t-"-mode and n_ particles in the "—"-mode, and a„_^_„_ 

are the respective amplitudes. For analysis of the SPDM it is convenient to introduce an alternative representation 
of the state in terms of the creation operators^ 

\4,)=Y,Wk,aUlm, (17) 

k,q 

where |0) is vacuum and Wkq is a symmetric matrix. Comparing Eqs. (|16p and (jl7p one finds w++ = ck2,o/'\/2, 
W-- = ao^2/V^ and — ai,i/2. Using Eq. ([T7|) one obtains the SPDM 

K^Aw^'w. (18) 

The SPDM is a 2 x 2 matrix and, therefore, its eigenvalues are completely determined by the determinant of the 
SPDM det[i^] — A1A2 and its trace Tr[i^] = Ai -|- A2 = with N being the number of particles. The entanglement, 
in turn, is determined by the normalized eigenvalues Ai^2 = which are found as 



Ai,: 




2C 



N 



(19) 



where |Cp = det K is the concurrence^ In order to see the relation with the standard definition of the concurrence 
in the two-particle case we introduce the "spin fiip" transformation ay |+) = — i |— ), ay |— ) — and the spin flip 

V'^ = {ay (g) ay lip))*, then 

C = = 4det[w]. (20) 



state 
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Thus the two-particle case completely fits into the canonical quantum-mechanical description. 
As follows from Eq. (I^Hl) the concurrence is expressed in terms of the amplitudes akq as 

C = 2a2,oao,2 - (21) 

For a two-particle state to be disentangled the amplitudes have to meet the condition C = 0. Let initially the state 
be disentangled. It is seen that this condition not necessarily holds for all t > if the dynamics of the system 
inhomogeneously depends on the population numbers. Since the Hamiltonian (jl5l) conserves the number of particles 
in each mode, the states |n+,n_) are the eigenstates and the time dependence of the amplitudes can be easily found 



_n_{t) = exp [-it (n+e+ + rt-f- + Ae„_|_,„_)] a„^,„_(0), (22) 



where Ae„^^„_ — [/+n+(n+ — 1) + U-n-{n^ ^ 1) + Substituting Eq. into Eq. (PT|) we find for the 

initially disentangled state 

C{t) = 2m2 i(o)e-»"* sm{ut), (23) 
where 17 = 2e+ + 2e_ + [/+ + [/_+ U+- and 

= [/_. (24) 

Several interesting conclusions can be drawn from this result. First, the entanglement oscillates between and the 
maximum value determined by the contribution of |1, 1) into the initial state. The origin of the oscillations can be 
traced to the structure of the concurrence, Eq. (PT|) . and the frequencies of the many-body amplitudes in Eq. (|22p . 
The interaction leads to the energy shifts Ae„m, which depends on the population of the particular modes. The phase 
mismatch between the amplitudes results in the nontrivial time dependence of |C(f)p. 

Second, the interplay between the effects of the intra-mode, cx Uu, and inter-mode, , interactions is not 

straightforward. Let the interaction between the modes be absent, = 0. As follows from Eq. ([25]) . even in this 

case initially disentangled states become entangled. 

The related effect is the mutual cancelation of the phase desynchronization if = U+ + U-, when despite the 

interaction, which changes the energies of the many-body states comparing to multiples of the single-particle states, 
initially disentangled states remain disentangled. 

The two-particle case is useful for establishing the relation with the standard description of entanglement. However, 
in order to grasp the general structure of the dynamics of entanglement of the quantum field with interaction it is 
constructive to consider more general case with an arbitrary number of particles N. In this case it is convenient to 
introduce Schwinger's angular momentum modeljii^ Defining the components of the operator of the angular momentum 
in terms of the creation and annihilation operators as 

Jx = ^(a+a- + 0-0+), 

Jy = ^(«+«- - 0-0+)> (25) 
Jz = ^(a+a+ - ala-), 

one can enumerate the states by the total angular momentum j — (n_|_ -t- n^)/2 and its projection m = (n+ — n^)/2 
instead of the population numbers n+ and n_, that is — \ j,m)g. 

Expressing the SPDM in terms of the mean values of the operator of the angular momentum one finds the concur- 
rence C{t) = N"^ /A — I (J(i)) P, so that the normalized eigenvalues of the SPDM can be expressed as 

\.2^\{l^J), (26) 
where J — 2\ (J) |/A^. Thus, the entanglement can be written as En — F{J), where 

^(^) = ^ E [1 + (-1)"^] l°g2 [1 + (-1)"^] • (27) 

n=l,2 



In particular, disentangled states are the states with the maximum magnitude of the average angular momentum 
I (J) I = j and completely entangled ones are those with | (J) | = 0. 
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In terms of the operator of angular momentum Hamiltonian (|T5|) can be presented as 

n = h + hJz + ^J"^. (28) 

where /i ^ MY.k ekUk/2 + Af^{U+ + U- + C/+-)/4 and /a = e+ - £- + {U+ - [/_)(! + Af) depend on the operator 
of the total number of particles, Af ~ a^a+ + oLa-, and are irrelevant for the entanglement dynamics and w is given 
by Eq. (IMI)- 

As follows from Eqs. ([26|) and psp the dynamics of entanglement is subject to the general constraint {Jz{t)) — 
{>Jz{0)), which follows from [H, Jz] — 0. Thus, the variation of entanglement is determined by the change with time 
of the "transversal" component J2(t) = {Jx{t)f + {Jy{t)f = (J+it)) {J-{t)), where J±{t) = J^{t) ±ijy{t). In 
particular, this imposes the upper limit on the value of entanglement produced by the interaction, Eisi{t) < E"^"^-* = 

F[{jzm/3]- 

Let us consider the dynamics of initially disentangled states. As follows from Eq. (pS)) . any such state can be 
obtained by rotating the state 

\i^[Pi,fi2.h)) = exp(-zJ,/3i)exp(~zJ'j,/32)exp(-zJ,/33) 1^,^)^, (29) 

where Pi are the Euler angles. Indeed, the average (J) transforms under rotations as a 3D vector. Hence, by rotating 
around z- and y-axes we can set in the new frame only J7z ^ 0. Since disentangled states have the maximum magnitude 
of the averaged angular momentum, this implies that in the new frame the state is proportional to 

Rotations around the z-axis do not affect and, therefore, the dynamics of entanglement of the states related 
through such rotations are identical. Hence, it suffices to consider only the states |'0(/3)) — |V'(0, /3, 0)). From 
[H; Jz] = it immediately follows that 

{Jz{t))p=jcos[li), (30) 

where the index (5 emphasizes the structure of the initial state, i.e. (. . ) ^ = ■ ■ ■ ■ The nontrivial part of 

the time dependence of {J±{t)) ^ is given by the last term in Eq. (|28t because the first two terms yield only the phase 
factor, which does not contribute to J^{t). Next, we notice that the solutions of the operator equations of motion 
j±{t) = i[Ho,J±{t)] with Hq = ujj^ have the form 

J± {t) = exp {-iijt ± 2iujJzt) J± . (31) 

Using these solutions we obtain (see Appendix |X] for the details) 

P =cos^(/3)+ 

+ 4sin2(^/2)cos2^^2(g,/2) (32) 
X |e''^sin(;3/2)sin(6l/2) + cos(/3/2)cos(6i/2)|^ , 

where sin(6'/2) = sin(/3) sin(a;i) and cot(7) = cos(/3) tan(a;t). The overall dependence of J on time and on the structure 
of the initial state parametrized by the angle p is shown in Fig. [T] 

It follows from Eq. that for /3 ^ or tt the transverse component of the effective angular momentum oscillates 
with the period T = 2t: /u. The maxima of J, which, according to Eq. (j26p . correspond to the minima of entanglement, 
are reached at = Tn/2 with integer n. As follows from Eq. (j32p . J(tmin) = 1 yielding Ejf — in agreement with 
the definition of the states |'0(/3)). There are also two sets of maxima of the entanglement at tmax = (1/4 + n)T and 
t'^'ax = (3/4 + n)T, where one has 

J\€L) = cos2(/3) +4sin2(/3/2)cos2(3/3/2)cos2^-2(^)^ 

respectively. One can see that J^(imaa;)/J^(ima2;) > 1. With increasing the number of particles the ratio monotonously 
tends to 1. 

In order to better understand the effect of the number of particles on the time evolution of entanglement we consider 
the case (5 — 7r/2, that is the initially disentangled states with symmetrically populated modes ((V'l Jz = 0). As 
follows from Eq. ([5^ and (|26l) these are the only states (among all initially disentangled states) that yield maximum 
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FIG. 1. Dependence of J{/3,t) (the concurrence) of initially disentangled two-particle states of the two-mode self-interacting 
boson field on time and on the structure of the initial state. 



entanglement Ej^ = 1 in the course of the time evolution. It can be shown that all such states are spanned by two 
orthogonal states corresponding to /3i = and 7r/2 in Eq. ((29|) . From Eq. (|32|) one finds that for these states 

j2 = cos2^-2(,^t). (34) 

The time dependence of entanglement following from Eq. (p4l) has two specific features. First, this is the periodic 
function of time. The entanglement evolves from 1 to and back within the vicinities tmin- When the number of 
particles increases these regions narrow. It is constructive to consider the limiting form of the time dependence when 
the number of particles becomes very large. In the limit iV 3> 1 one can approximate 

oo 

j2(i)w^exp[-(A^-l)(t^t-7rn)2]. (35) 

Thus, away from the points where J(t) = 0, it can be regarded as the sequence of the Gaussian bumps localized near 
tmin- As a result, when ^/N ^ 1 one can consider the system as being in the completely entangled state most of the 
time. ^ 

The Gaussian decay of J{t) is similar to the Gaussian decay of the coherence of the central system,^-' two spins 
1/2 coupled to the bath. In particular, the same dependence of the decay rate on the number of particles cx ^/N 
in the environment [notice — 1 in Eq. (I35p ] should be emphasized. There are, however, two important differences 
between this situation and our case. First, the Gaussian decay of the coherence of the central system appears in 
the limit of the slow dynamics of the environment. In the opposite limit the decay follows the Lorentsian law and 
in the intermediate case both types of decays present at different time scalesJ^ For the two-mode boson field it is 
meaningless to separate particular particles and environment due to indistinguishability, however, it is worth noting 
that J{t) does not depend on the single-particle energies. The second important difference is that the concurrence of 
the two-mode boson field exhibits oscillations while the loss of coherence of the central system is irreversible. 



IV. ENTANGLEMENT DYNAMICS IN THE JAYNES-CUMMINGS MODEL 



In the previous section we considered two basic features of the entanglement dynamics in the system of coupled 
quantum fields — redistribution of the initial total entanglement and the production of entanglement by the interac- 
tion. In the present section we apply these ideas for the analysis of entanglement in a more complex situation. 

We consider a single two-level atom interacting with the quantized electromagnetic field. The electron transitions 
in the atom are assumed to be characterized by definite helicity. The excitation of the electron state with the spin 
down at the ground level, \g i), into the spin up state at the excited level |e f) occurs through the absorption of 
-polarized photon and so on. The dynamics of the system is described by the Hamiltonian of the two-mode 
Jaynes-Cummings (JC) model^ which we write down in terms of the creation and annihilation operators 

n = J2 4^^4c« + E 4''4«fc + n+ + n^. (ae) 

K k 
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Here the first two terms describe the dynaniics of the free atom and the free field, respectively, with k and k running 
over the atomic states, {g l,e f, e i}, and the photon polarizations, + and — , respectively. The interaction 
between the atom and photons is described by the Hamiltonians Hi = uJiaiO-j + h.c, where cr^_ = cl^Cg^, al^ = cl^Cgf 
and u!± are the respective Rabi frequencies. 

As the consideration in the previous sections suggests, the main effect on the entanglement dynamics is due to 
the coupling between the systems rather than due to their internal dynamics. However, the interplay between 
different characteristic frequencies may lead to appearance of beatings. Therefore, in order to eliminate the technical 
complications and to concentrate on the main features, we adopt the resonant approximation and set = e)f = 0. 
Additionally we assume that the symmetry between transitions with different helicities is not broken so that a;+ = 
a;_ = UJR. 

We consider the photon entanglement, which is defined as the von Neumann entropy of the single-photon density 
matrix with the matrix elements 



k,{t) = (4W««W) • (37) 



This is 2 by 2 matrix and the problem of entanglement can be approached using the same description as in the 
previous section. As well as before our main objective is to study the time evolution of entanglement of initially 
disentangled photonic states. This closely corresponds to the situation when, for example, the cavity in the ground 
state is pumped by an external source. 

Taking into account that initially the atom is not excited the state of the system can be presented as 

|*(0)) = |V(/3))|0),, (38) 

where is the disentangled photon state obtained by rotations from the state |j, [see Eq. (|29|) ] and |0)g 

denotes the state of the atom with both electrons at the ground level. 

The time dependence of the photon entanglement of the states with only a few photons is highly nontrivial comparing 
to the few particle case of the two-mode boson field considered in the previous section. Here we would like to note only 
one characteristic feature leaving more detailed discussion for a separate publication. Drawing the analogy with the 
consideration of the two-particle case in the previous section one can expect that in the present case as well the time 
dependence of the entanglement will be determined by the mismatch between the amplitudes of states with different 
population numbers. In fact, as will be evident shortly, the typical frequencies are determined by the square roots 
of the population numbers. Taking this circumstance into consideration one can expect that the time dependence 
of the concurrence is the result of superposition of several harmonics with incommensurate frequencies. Thus it is a 
quasi-periodic function with a complex profile. 

With increasing the number of photons, however, the contribution of specific frequencies becomes less important 
and the overall shape of J{t) changes toward some general regular pattern shown in Fig. [21 In order to describe it 
we use the explicit form of the Heisenberg representation for the photon operators ak{t) — exp(ii'Hfc)afc exp(i<'Hfc). 
Taking into account the separability of the dynamics we have^ 



ak{t) = e 



1 [Cktj - iCk- 



sin (Ckt) 1 , , 

ak - 1 y ^yk , (39) 



where = wjjo^afc and = U!ii{a\,ak + 1). It should be noted that the second term in Eq. (1391) does not contribute 
to Kkq{t) because it vanishes while acting on the initial state of the atom. 

First, we consider {Jz{t))- Taking into account the identity / (C^) ak — akf {Ck), which holds for any well-behaving 
function /, we find 

iJ.it)) = iJ.) + ([cos(2C_i) - cos(2C+i)]) , (40) 

where the average is taken over the initial photons state, (...) = (V'(/?)l ■ • • IV'(/?))- The second term in the r.h.s. of 
Eq. (Uni) plays significant role only in the case when all photons have the same polarization, that is when the effect 
of the interaction is obviously small. In this case the interaction produces small entanglement ^ 1 /N oscillating with 
the frequency Nujfi,. 

In the unpolarized case, when (J7z) = the last term in Eq. (|40p vanishes identically because of the symmetry of 
such disentangled states with respect to flipping all photon spins. While in the weakly polarized case the ratio of the 
last term to (J^z) can be shown to remain limited from above and to decrease with the number of particles oc 1/\/N. 

Thus, while due to absorption of the photons {Jz{t)) is not a constant, as it was for the two-mode boson field, 
its total variation is small in the limit iV ^ 1. Therefore, the initial value {J'z{0)) can be considered as imposing a 
limitation on the highest entanglement, which can be reached. 
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FIG. 2. Time dependence of J{t) for the Jaynes-Cummings model. The time is measured in inversed coupling constant ujh. 
(a) Comparison of the exact dynamics (dotted lines) with the approximation defined by Eq. (|42|l (dashed lines) on the long 
time scale for the system with 50 (lower panel) and 200 (upper panel) photons, (b) The short time dependence for 10 particles 
— the exact result (solid line) and the approximation (dashed line). The transition to the Gaussian law occurs fairly quickly 
with the number of particles. It is illustrated in inset, where the short time dependence is shown for 200 photons in log-scale. 



Because of these reasons, we consider m more details the case when (J7z(0)) = 0, that is when the initial photon state 
is given by Eq. ([29|) with /32 = tt/2. In this case the eigenvalues of the photon SPDM, and, hence, the entanglement, 
are determined by the magnitude of (j7+(t)). In order to calculate it we assume that the main contribution is due to 
the states with small m, that is n+,n_ ^ 1 and n+ w n_. This approximation improves with increasing the number 
of particles. 

The best way to employ this approximation is to derive the effective photon Hamiltonian directly from Eq. p9p . 
Taking into account that initially the atom is in the ground state we can effectively substitute Eq. p9p by ak(t) = 
exTp{iCkt)ak exp(— iC^t), which is the Heisenberg representation induced by the Hamiltonian 

H = L^R^j + Jz + lob.Vj-Jz. (41) 

Using the assumption regarding the main contribution due to small m we can expand this Hamiltonian with respect 
to J^z/j finding in the first nonvanishing order 

^--«V7-^^^ (42) 

Comparing this expression with Eq. (|28|) we come to the conclusion that in this limit the photons behave as the 
two-mode self- interacting boson field with the effective interaction strength w = lor/^.j'^/'^. The characteristic feature 
of the effective interaction is that its intensity decreases with the number of particles oc l/N^/"^. This "spreading" 
is caused by sharing the interaction with the single atom among all photons and leads to different dynamics of J{t) 
comparing to the one we ha ve see n in the previous section. Using the results obtained there we find the period of the 
long-scale oscillations T — \'2N^ /lor and the decay time r — V^N/ujR (see Fig. Hi). 

It should be noted that, as can be seen from Fig. [51 the shape of the bumps changes with their number — they become 
asymmetrical and acquire the side oscillations, the revival of coherence occurs non-monotonously. The approximation 
we have used misses these long-time changes, which are the reminiscent features of the quasi-periodicity mentioned 
above. If, however, one restricts to the initial growth of entanglement (the drop of coherence) these features are not 
important, and one has in the limit A'^ ^ 1 

J(i) = e-("«*/2W)'. (43) 



V. CONCLUSION 



We have considered the basic dynamics of entanglement in the system of coupled quantum fields. We have quantified 
the entanglement as the von Neumann entropy of the single particle density matrix (SPDM) and studied its time 
evolution for different situations. 

The first question, which has to be answered is how is it possible to produce an entangled state of a quantum 
field. The circumstance, which makes this question non-trivial, is that by classical source entangled states cannot 



10 



be excited. Following the straightforward analogy with the canonical quantum mechanical picture one can consider 
the system of coupled quantum fields. We show, however, that if the coupling is linear, i.e. if the total Hamiltonian 
is a quadratic function of the field operators, the total entanglement conserves and one, at most, can only have the 
transfer of the entanglement between the coupled fields. 

The simplest system demonstrating entanglement of initially disentangled states is the two-mode boson field with 
self-interaction. Reformulating the problem using the formalism of Schwinger's model of angular momentum we show 
that there is the direct relation between J, the magnitude of the average angular momentum (J) and the concurrence. 
More specifically, the states with maximum possible magnitude of the average angular momentum are disentangled 
and those with J = are completely entangled, meaning that the SPDM is proportional to the identity matrix. 

We show that in the limit of high number of particles J{t) has the form of the periodic sequence of Gaussian bumps, 
whose width (inverse entanglement time) decreases with the number of particles oc 1/^/N. The interesting feature 
is that the entanglement is produced even if the interaction is between the particles within the same mode but not 
between different modes. This phenomenon is in the striking contrast with the standard quantum mechanical picture, 
where separable dynamics cannot entangle initially disentangled particles. The origin of this effect lies in the structure 
of the many-body states of the field. Because of indistinguishability of the particles one cannot say which particle 
belongs to which mode. As a result only states with amplitudes meeting the special condition are disentangled. The 
dynamics nonlinearly depending on the number of particles may break this condition resulting in entanglement. 

We apply these results for analysis of the two-mode Jaynes-Cummings model. The photon-atom interaction has 
been assumed to be helicity preserving thus the dynamics of "-|-"- and "—"-polarized photons are separable. However, 
absorption and re-emission of photons by the atom introduces an effective interaction between the photons of the same 
polarizations. This interaction, similarly to the case of the two-mode boson field, entangles photons with the typical 
Gaussian drop of coherence in the limit of large number of photons ^ 1. Since the interaction with the single atom 
is shared among the photons the strength of the effective interaction drops cx l/N^^^. This leads to prolonged both 
the oscillations of entanglement T = \/2N^/u!ii and the entanglement time t = . 
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Appendix A: Dynamics of the transversal angular momentum 



For the different situations considered in the main text we need the value of the typical matrix element 

J+(0,/3) = ii^m eM^Jz<f>)J+ , (Al) 



where |V'(/3)) = exp{-ijy^) \j,j)g. 
Using the transformation rule 



we can rewrite Eq. (jAip as 



^■^JyPj^e-^JyP = sin(/?) + J+ cos2(/3/2) - sin2(/3/2). 



J+(0,/3) = s {j,jmJz sin(/3) - sin2(/3/2)] |j, j)^ , 



(A2) 



(A3) 



where we have taken into account that j)5 = and have introduced TZ = exp(ij7j,/3) exp(ij'z(l>) exp(— zj7y/3) the 

operator of rotation by angle (j) around the axis along the direction (sin(/3), 0, cos(/3)). Expanding TZ in terms of the 
Euler angles (7', 6*, 7) related to (3 and through 



sin(6l/2) ^ sin(/3)sin(<?!)/2), 
tan[(7' + 7)/2] = cos(^) tan(0/2). 



7-7 = TT, 



we find 



J+(0,/3) = e''('^'+'') \jsm{/3)d,,,{0) - e^^V^ sin2(/3/2)d,_i,,(0). 



(A4) 



(A5) 
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Here 



dm AO) 



(j + TO)!(j - m)! 



cos^'+™(0/2)sin^-"(e'/2) 



(A6) 



are the respective matrix elements of the {2j + l)-dimensional irreducible representation of exp{—ij^y9)}^ Using 
Eq. (M\i in Eq. (K5\ we finaUy obtain 

J+{(j), 13) = 2je-*('''+'^' sin(/3/2) cos2^'-i(6i/2) [cos(/?/2) cos(6i/2) - e''' sin(/3/2) sin(6i/2)] . (A7) 

In the case of the special interest /3 — tt/2 this expression significantly simplifies 

J+(</.,7r/2) = je-^*/2cos2^"-i(0/2). (A8) 

For small (j) and large j the magnitude of this function decays with 4> following the Gaussian law. Let us denote 
J{t) = cos^J-\(j)/2), then 



^ = i(2j-l)tan(0/2)JW. 
For ^ 1 we can keep only the leading term in the Taylor expansion of tan((/)/2) and obtain 

1 ^-^^ 



J(i)«exp --(2j-l)(^| 
Such representation is meaningful for sufficiently large j when j(fp > 1 



(A9) 
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